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Abstract 

Models of neutrino mass which attempt to describe the observed lepton mixing pat- 
tern arc typically based on discrete family symmetries with a non-Abclian and one 
or more Abelian factors. The latter so-called shaping symmetries are imposed in 
order to yield a realistic phenomenology by forbidding unwanted operators. Here we 
propose a supersymmetric model of neutrino flavor which is based on the group T-j 
and does not require extra Z^v or U(l) factors, which makes it the smallest realistic 
family symmetry that has been considered so far. At leading order, the model pre- 
dicts tribimaximal mixing which arises completely accidentally from a combination 
of the T-j Clebsch-Gordan coefficients and suitable flavon alignments. Next-to-leading 
order (NLO) operators break the simple tribimaximal structure and render the model 
compatible with the recent results of the Daya Bay and Reno collaborations which 
have measured a reactor angle of around 9° . Problematic NLO deviations of the other 
two mixing angles can be controlled in an ultraviolet completion of the model. 
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1. Introduction 



The triplication of chiral families remains one of the biggest mysteries in particle physics. 
A clue towards unraveling the principle behind the fact that quarks and leptons come 
in three copies is provided by the observation of a remarkable lepton mixing pattern: 
contrary to the quark sector, the mixing of the leptons is described by two large and 
one small angle. Until recently, the results of neutrino oscillation experiments were well 
compatible with a PMNS matrix of the intriguingly simple structure [1-3] 

/ 2 1 n \ 

(1) 

This so-called tribimaximal mixing suggests an intimate connection of the three genera- 
tions of leptons, which can be realized in the framework of non-Abelian discrete family 
symmetries. Imposing such a horizontal symmetry G allows to unify different genera- 
tions into a multiplet of the given non-Abelian group. With three families, the physically 
interesting groups should have a triplet representation, limiting possible choices to sub- 
groups of U{3). Many of these have been successfully applied to construct models of 
tribimaximal lepton mixing, see for instance [4-6] and references therein. 
In general, family symmetry models can be classified according to the origin of the 
symmetry of the neutrino mass matrix. We will assume that neutrinos are Majorana 
particles; then their mass matrix is always symmetric under a Klein symmetry Z2 x Z2. 
Working in a basis where the charged leptons are diagonal,^ the explicit form of the 
Klein symmetry expressed in terms of 3 x 3 matrices is dictated by the PMNS mixing 
matrix and can be determined as 
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Here p = 0, 1 and g = 0, 1, yielding a symmetry group of four elements. This neutrino 
flavor symmetry can arise as a residual symmetry of the underlying family symmetry G, 
in other words, the four elements of Eq. (2) can also be elements of the imposed family 
symmetry. Models of this type are called direct models [14]. In indirect models, on 
the other hand, the above Klein symmetry is not a subgroup of G. Models of this 
class are typically based on the type I seesaw mechanism [15-18] with the assumption 
of sequential dominance [19-21], i.e. a hierarchy among the three terms arising from 
three right-handed neutrinos. Here, the main role of the family symmetry consists 
in explaining special vacuum configurations of the flavon fields that break the family 
symmetry [14]. In addition to these pure classes of models there are semi-direct models 

^This classification can also be applied to GUT models where the charged lepton mass matrix is 
typically only approximately diagonal. In such a setup, the total PMNS mixing matrix will involve 
charged lepton corrections that have to be taken into account separately, see e.g. [7-11], leading to 
characteristic mixing sum rules [12, 13]. 
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in which one of the Z2 factors of the Klein symmetry arises as a residual symmetry of G, 
while the other factor arises accidentally. In fact, the famous Altarelli-Feruglio model 
based on A4 [22, 23] belongs to the semi-direct class as the tribimaximal /i-r symmetry 
is not part of A4. 

Hitherto, regardless of the type of model, the non-Abelian family symmetry has always 
been augmented by extra Abelian factors such as Z^r and U(l) in order to yield realistic 
phenomenology. These shaping symmetries were crucial for controlling the coupling of 
the Standard Model neutral flavon fields to the leptons. In [22,23], for instance, the 
neutrino and charged lepton sectors are separated by means of a Z3 shaping symmetry. 
In that sense, one should speak more precisely of the Altarelli-Feruglio A4 x Z3 model, 
and the family symmetry G is thus a group of order 36, rather than 12. 
Mindful of this subtlety of defining the full family symmetry of a model, a systematic 
scan over 76 different groups has been performed in [24] with the purpose of studying 
whether or not there is an inherent connection between ^4 and tribimaximal mixing. The 
results of the scan proved that there are indeed several, even more minimal, groups which 
are capable of describing tribimaximal mixing (assuming the imposed simple alignment 
of the flavon vevs can be justified). The smallest such group was identified to be T7, a 
group of order 21 which is sometimes also called the Frobenius group Z7 x Z3 [25-27]. 
Note that this non-Abelian symmetry does not require any extra shaping symmetry. 
It is furthermore interesting to point out that T7 has no Z2 subgroups, and hence the 
Klein symmetry of the neutrino mass matrix cannot be part of this family symmetry. 
As such, a corresponding T7 model would be neither direct nor semi-direct, but rather 
of indirect type. Not necessarily requiring the seesaw mechanism, it would - to the best 
of our knowledge - be the first indirect model which is not based on the assumption 
of sequential dominance and the quadratic appearance of flavon fields in the effective 
neutrino mass term. 

It is the purpose of this article to present the details of a complete T7 model of leptons 
which yields tribimaximal mixing at leading order, including a discussion of the vacuum 
alignment. We put particular emphasis on explaining how the Klein symmetry arises 
completely accidentally from a combination of the T7 Clebsch-Gordan (CG) coefficients 
and suitable flavon alignments. As tribimaximal mixing has been ruled out by the 
Daya Bay [28] and Reno [29] measurements of a reactor mixing angle ^13 of around 
9°, it is necessary to investigate if next-to-leading order (NLO) effects can generate 
large enough deviations from the tribimaximal leading order prediction. With the Klein 
symmetry arising purely accidentally, NLO terms are bound to perturb the tribimaximal 
mixing pattern. While switching on ^13, NLO corrections will in general also give rise 
to perturbations of the atmospheric and, more critically, the solar mixing angle which 
may be phenomenologically unacceptable. However, these unwanted NLO corrections 
can be suppressed in an ultraviolet completed version of the T7 model involving special 
messenger fields. 

The remainder of this article is organized as follows. We define and discuss the minimal 
Tj model in Section 2, both at leading as well as next-to-leading order. The ultraviolet 
completion of the model is presented in Section 3. Section 4 addresses the question 
of vacuum alignment, and the conclusions are drawn in Section 5. The relevant group 
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Table 1: The transformation properties of the minimal Ty model of leptons. 

theoretic details of T7 are laid out in Appendix A, including the generators of the rep- 
resentations in a basis with a diagonal order three element as well as the corresponding 
CG coefficients. 

2. The minimal T7 model 

2.1. Tribimaximal mixing at leading order 

In this section we describe a minimal model of leptons based on the family symmetry T7. 
This non-Abelian finite group comprises 21 elements and has three singlet as well as two 
triplet representations, which we denote by 1, 1', 1", 3 and 3, respectively. As the 
model does not feature any shaping symmetries, its structure is solely determined by 
the family symmetry T7 as well as the gauge symmetry SU(2)2, x U(l)y. We work in a 
supersymmetric framework with two Higgs doublets, hu and /i^, transforming trivially 
under T7. The three generations of left-handed lepton doublets L are unified in a triplet 
representation of T7, while the right-handed charged leptons e^, fi^ and live in the 
three distinct one-dimensional representations. In order to break the family symmetry, 
two fiavon fields are introduced, namely if and transforming as a 3 and 3, respectively. 
The particle content of this model is summarized in Tab. 1. For later purposes we have 
also listed the hypercharges as well as the charges under the U{1)r symmetry which is 
required for the F-term fiavon alignment mechanism discussed in Section 4. 
With these T7 assignments, the leading order superpotentials of the charged lepton and 
the neutrino sector read 

W, = y^fLe^h. + y^lLfi^h. + yrfLr'^h,, (3) 
Wu = yi^LLhuhu + y2j2LLhuK , (4) 

respectively. Here, the y's are dimensionless coupling constants, and A denotes a common 
cut-off above the scale of family symmetry breaking. The neutrino mass terms originate 
from the Weinberg operator [30] augmented by a flavon field. 

It is important to note that only one fiavon, ip, enters the charged lepton sector at leading 
order. The other fiavon, ip, cannot couple to the charged leptons due to the absence of a 
singlet representation in the T7 Kronecker product 3 ® 3, see Appendix A. We will see 
below that the vev of ^ can be chosen such that it breaks T7 down to Z3, thus leading to 
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a diagonal charged Icpton mass matrix. In the neutrino sector, both flavons are present, 
as the symmetric product of 3 ® 3 contains both a 3 and a 3. 

In order to find the vacuum alignments that lead to tribimaximal mixing, we insert the 
flavon vevs in Eq. (3) and Eq. (4) in their most general form. Contracting the Ty indices 
using the CG coefficients given in Appendix A and inserting the Higgs vevs Vu and v^, we 
obtain the respective mass matrices. In the left-right convention for the charged leptons, 
where the left-handed particles are to the left of the Yukawa matrix, we find 
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This simple pattern suggests the following fiavon vacuum configuration to obtain tribi- 
maximal mixing: 



'V 
.0> 



(7) 



The alignment of (f ensures a diagonal charged lepton mass matrix in Eq. (5), and brings 
the first line in Eq. (6) into tribimaximal form. The alignment of (p is of an extremely 
simple form and generates a contribution to M^, which, again, has tribimaximal structure. 
Adopting the flavon alignments of Eq. (7), our minimal T7 model therefore predicts the 
following charged lepton and neutrino mass matrices at leading order. 
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implying a tribimaximal PMNS mixing matrix and the complex valued neutrino masses 
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Before continuing to the study of NLO effects, let us pause for a moment and appreciate 
the beauty of the leading order model. With regard to the charged lepton sector, we 



4 



observe that the ahgnment of (p is left invariant by the diagonal order-three generator 
d of T7, see Tab. 5 in Appendix A. The corresponding Z3 symmetry which is thus 
preserved in the charged lepton sector after the family symmetry is broken enforces 
a diagonal M^. In contrast, the tribimaximal Klein symmetry Z2 x Z2 of My arises 
completely accidentally from a combination of the T7 CG coefficients and suitable fiavon 
alignments (which will be derived and thus justified in Section 4). This situation is 
reminiscent of the Altarelli-Feruglio model [22,23] where one Z2 is contained in Aa, while 
the other Z2, corresponding to a /i-r symmetry, arises accidentally. In the case of T7, none 
of the two Z2 factors of the tribimaximal Klein symmetry is a part of the underlying 
family symmetry. As a consequence, one expects no protection of the tribimaximal 
structure when NLO terms are taken into account. 



2.2. Next-to-leading order effects 

The next-to-leading order superpotential for the charged and neutral leptons contains 
all terms up to mass dimension six and seven, respectively, that are invariant under the 
gauge and family symmetries. Here, we only indicate the terms that are not already 
part of the leading-order superpotential: 
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(12) 



The subscripts in Eq. (12) indicate that we have contracted the family indices of the 
corresponding product to obtain a 3 and 3, respectively. As before, contracting the 
family and gauge indices, and substituting the vevs for the Higgs and fiavon fields gives: 
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A quick inspection of Eq. (13) shows that the contributions proportional to Cf , Cf , CI 
are identically zero and have no effect whatsoever on the mass matrices. The contribu- 
tions proportional to C|, C3 , do not disturb the diagonal structure of Mi and hence 
do not affect the bi-unitary transformations that diagonalize it and that enter in the 
definition of f/pMNS- In contrast, Cg, , will lead to a departure from tribimax- 
imal mixing, where it has to be noted that the magnitudes of the coefficients have to 
be weighted by the corresponding Yukawa couplings, and hence the effect of C| will be 
neghgible. 

From Eq. (14) we see that the term proportional to Ci is identically zero and will not 
have an impact on the mixing matrix. The term corresponding to C3 is non-vanishing, 
but commutes with the two generators of the tribimaximal Klein symmetry (see e.g. [4] 
for the explicit form of the generators), and as a consequence does not disturb the 
tribimaximal form of the mixing matrix. 

The NLO terms that do cause a departure from tribimaximal mixing are those corre- 
sponding to 6*2, C4, C5, Cg. The contribution proportional to C4 is invariant under 
a /i — r symmetry (which corresponds to one of the generators of the tribimaximal 
Klein symmetry) and thus leaves 6*23 and 6*13 unchanged. Of particular interest is the 
term proportional to Cg , since its effect is such that 6*12 is almost unchanged, 6*23 stays 
within its 3a interval, and 613 receives a sizable contribution. Therefore it seems to be 
a good strategy to try and suppress all NLO contributions to the superpotential except 
the one proportional to C5 . Thus, we will keep the deviation of 6*12 and 6*23 from the 
tribimaximal case small, and at the same time achieve a non-vanishing ^13 in agreement 
with experiment. In Section 3 we present a renormalizable, ultraviolet completion of 
our model that meets exactly these criteria outlined above. Before that, however, we 
will explore in Section 2.3 the phenomenology of our model at the effective level. We 
first consider the model in its most general form without any specific assumptions about 
the coefficients. At the end of Section 2.3 we also discuss the case where all coefficients 
which break the tribimaximal structure, except for C5 , are set to zero. 

2.3. Phenomenology 

We now compare our results with experiment [28,29,31-33]. To that end, we first 
have to specify the coefficients in the leading and next-to-leading order superpotential 
(cf. Eq. (3), Eq. (4), Eq. (11), Eq. (12)). The mixing angles, CP phases and masses will 
then be uniquely determined. 

The leading-order superpotential gives tribimaximal mixing which is form-diagonalizable 
[34]. As such, the mixing angles are in principle independent of the coefficients i/e, y^, 
Vt, yi, 2/2- However, changing the coefficients affects the neutrino masses and can lead 
to a reordering of the mass eigenstates that in turn has an impact on the mixing angles. 
Since we work in a basis where the charged lepton mass matrix is diagonal, we can 
identify ye, y^, yr with the corresponding Yukawa couplings for which we substitute 
their experimentally determined values [35]. For yi and y2, we choose values such that 
we can fit the mass ratio Am'^i/Am'^i. 
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Figure 1: Contour lines of Am^j^/Aml^ = 30 as a function of yiv^ and y2V^p for normal neutrino 
mass ordering. The green region corresponds to values of the parameters where the mixing angles are 
tribimaximal. 

In Fig. 1 we present the contour lines of Am|^/Am2i = 30 as a function of yiv^p and 
y2Vip (cf. Eq. (9) on page 4) for normal neutrino mass orderings. Here we do not consider 
the case of inverted neutrino mass orderings as it would require yif^ <^ y2V,p. The green 
region in Fig. 1 corresponds to values of the parameters for which we obtain tribimaximal 
mixing, whereas in the red region, the mixing matrix corresponds to a permutation of 
the columns of f/pMNS? i-e. tribimaximal mixing and the desired mass hierarchy cannot be 
simultaneously satisfied. Note that we have chosen yi and y2 to be real for presentational 
purposes only. In the general case, yi and y2 are complex, and to fit the mass ratio is 
equally easy. 

In the next-to-leading order superpotential, we assume that the fundamental scale of 
the theory is about one order of magnitude larger than the flavon vevs and fix the ratio 
at v^/A = Vf/A = 0.10 (see Eq. (13) and Eq. (14)). For the coefficients C]", . . . , CJ, 
we choose 1000 sets of complex numbers whose real and imaginary parts are uniformly 
distributed between -2.5 and +2.5. For each such set, we calculate the corrections to 
the charged lepton and neutrino mass matrices and recalculate the mixing angles and 
masses. Note that now. Mi is not necessarily diagonal, and we have to resort to the 
general definition of ?7pmns in terms of the bi-unitary transformations that diagonalize 
the charged lepton and neutrino mass matrices. 

We present the results in Fig. 2(a) on the following page. The LO order prediction has 
given way to a "blob" centered around the tribimaximal values for the mixing angles. 
Remarkably, 99% of the points lie in the 3a interval [31] of 623, whereas 13% also satisfy 
the 3a bounds for 6*12 (see first panel of Fig. 2(a)). As was to be expected, the strongest 
constraint comes from ^13. From the third panel of Fig. 2(a) we see that low values 



7 



- 


















- 


- 


















- 
















































IP 














■ ? 




























































































10 20 ;! 


T 4 


n r,n no 70 so 9q 



- 














- 


















































i 


2 


1 30 


40 r, 


m 


^- 

n 7 




n m 



in 20 3n 40 r,o 00 to so so 
^23 



(a) t><^/A = v^/A = 0.10 and -2.5 < ^HeCf , 3m Cf < 2.5. 



10 20 ;io 40 go 70 so <)o 



10 20 30 40 .'jO go 7n so 90 




in 2n 3n 40 50 go to so 90 



(b) v^/A = 0.25, w^/A = 0.05, -2 < 5KeC^, amCg^ < 2, all other coefficients set to zero. 



Figure 2: The mixing angles for a given set of leading order coefficients yi, 7/2, Ve, Uii, Vt and 1000 sets 
of random values for the next-to-leading order coefficients Ci, . . . , CJ. Our choice for the coefficients 
Ui = 0.50, 2/2 = — 1.52 j/i Wi^/wip correspond to Arri^i/ Am^i ~ 30 in the LO approximation. The green 
bands correspond to the 3<j intervals [31] for the mixing angles (normal hierarchy in the case of 6*23 and 
^13). In each diagram, the color of the markers indicates whether the data point lies in the 3a interval 
of the respective angle that is not plotted (blue) or not (red). 



for 9i3 are preferred, but 10% of the points can accommodate a large reactor angle in 
agreement with experiment. This number drops to 1% when we take the constraints 
from the other two angles into account. 

In Fig. 2(b) we present the results for a model where all terms except that lead to 
departures from tribimaximal mixing are suppressed. This will correspond to the case 
of the renormalizable model that we will discuss in Section 3. As we can see from the 
first panel of Fig. 2(b), 612 is practically unchanged from its tribimaximal value, whereas 
623 varies very little, and both angles are well within their 3a error bands. At the same 
time, we can easily obtain a large ^13 in such a way that the experimental constraints 
on the other two angles are always satisfied (blue points in the third panel of Fig. 2(b)). 
At the same time we can observe an interesting correlation between 623 and 613, namely 
the larger ^13 is, the more 623 deviates from maximal mixing.^ Unfortunately, this effect 

^This numerical observation is a reflection of a well-known mixing sum rule relating the deviations of 
the atmospheric and reactor angles from their tribimaximal values via the CP phase S [36,37]. 
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is symmetric around 623 = 45° and hence does not give a hint whether this deviation is 
positive or negative. 

3. An ultraviolet completion of the model 

We have shown in Section 2 how the minimal T7 model, which predicts exact tribi- 
maximal lepton mixing at leading order, receives in general significant corrections at 
next-to-leading order. Such deviations from the tribimaximal pattern are welcome in 
order to accommodate the observed large value of the reactor angle 613. On the other 
hand, the same set of corrections can potentially modify the other two mixing angles 
to values disfavored by current global fits. In particular, the solar angle should obtain 
only small corrections, as its tribimaximal value of 612 ~ 35.3° fits the experimental 
result already remarkably well. In order to obtain a deviation from tribimaximal mixing 
which stabilizes the solar angle in first approximation, one can make use of corrections 
of the so-called trimaximal type [38-46]. These are typically defined by corrections to 
the tribimaximal neutrino mass matrix whose eigenvectors are proportional to either 
(1,1,1)"^ [37] or (2,-1,-1)"^ [47]. Comparison with Eq. (14) reveals that there exists 
exactly one NLO correction which respects the trimaximal structure but breaks the tribi- 
maximal one. This subleading contribution to the neutrino mass matrix is proportional 
to the coefficient and can be traced back to an NLO operator in Eq. (12) in which 
LL is contracted to a 3 of T7. Notice that the other two terms with LL contracted to a 
3, i.e. those proportional to Ci and C3 , do not lead to deviations from the tribimaximal 
structure. 

This begs the question of how to remove or suppress the unwanted NLO operators while 
keeping the one proportional to C5 . To this end, it is useful to recall that we have so far 
only discussed the model at the effective level, i.e. writing down all (non-renormalizable) 
terms which are allowed by the imposed symmetries. However, any particular ultravio- 
let completion of such an effective model does not give rise to all effective NLO terms. 
The underlying renormalizable model will involve messenger fields which have to be in- 
troduced to mediate the effective leading order operators. These necessary messengers 
typically do not generate all NLO terms [48]. If one particular NLO term is desired, as is 
the case in our Tj model, an extra messenger has to be introduced. In the following we 
adopt the strategy of considering an ultraviolet completion of the effective T7 model such 
that the operator proportional to C5 is allowed, while the bothersome NLO terms are 
forbidden or sufficiently suppressed. We begin the discussion of the ultraviolet comple- 
tion of the model with the charged lepton sector, before turning to the more interesting 
neutrino sector which will be responsible for the breaking of the tribimaximal to the 
trimaximal pattern. 

3.1. The charged lepton sector 

The leading order operators in Eq. (3) can be obtained by imposing a T7 triplet mes- 
senger pair B, with the transformation properties as given in Tab. 2. The resulting 
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Table 2: The charge assignments of the (effective) Higgs SU(2) ^ triplet A as well as of the messenger 
fields required for the ultraviolet completion of the minimal Tj model. 

renormalizable charged lepton superpotential takes the form 

Wp ~ LhdQ'' + Q^e" + e^/i" + Q^T^ + QQ^iMQ + + ^) , (15) 

where we have suppressed all dimensionless coupling coefficients. Note that 9 and 9'^ 
allow for a bilinear mass term. However, due to the absence of any shaping symmetry, we 
can also couple a flavon triplet to this product, thus leading to a trilinear term which, 
after family symmetry breaking, gives a correction to the messenger mass. This is a 
general feature for heavy messenger fields which transform as triplets under T7. Ignoring 
the correction to the messenger mass, we can integrate out the pair 9, 9^ and obtain 
the effective superpotential of Eq. (3) without any higher order corrections. Hence, 
in the present ultraviolet completion, the only source of deviations from a diagonal 
charged lepton mass matrix originates from the aforementioned trilinear terms correcting 
the messenger mass. Such a correction can be rendered sufficiently small by assuming 
<^ Mq. In the following, we will therefore ignore any NLO terms in the charged 
lepton sector. 

3.2. The neutrino sector 

We now turn to the neutrino sector where we want to formulate a renormalizable theory 
which, at leading order, gives rise to the effective superpotential in Eq. (4). The most 
popular possibility to derive the Weinberg operator from renormalizable terms is pro- 
vided by the famous type I seesaw mechanism [15-18]. However, applied to the minimal 
T7 model without any shaping symmetry, the type I seesaw mechanism does not yield 
the tribimaximal structure at leading order. This can be traced back to the fact that 
the triplet representations of T7 are complex. Putting the L into the 3 would entail to 
introduce right-handed neutrinos A^^ in the 3 so that the Dirac neutrino Yukawa term 
LhuN'^ leads to a Dirac mass matrix that is proportional to the identity matrix. Then 
the tribimaximal structure of the effective light neutrino mass matrix would be inherited 
directly from a tribimaximal right-handed neutrino mass matrix M/j. With N'^ being 
a 3 rather than a 3, the alignments of the 3 fiavon (f and the 3 fiavon (p would have 
to be exchanged to get Mji of tribimaximal form. However, changing the alignment of 
the fiavon coupling to the charged leptons to (1, 1, 1)"^ destroys the diagonal charged 
lepton mass matrix in Eq. (5). As a result, under the fairly general set of assumptions 
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Figure 3: The diagrams contributing to the leading order neutrino superpotcntial of Eq. (4). 



that we made above, the total leading order lepton mixing would no longer feature the 
tribimaximal pattern. 

As an alternative to the type I seesaw mechanism, we can adopt the type II seesaw [49- 
53]. For the sake of clarity, we assume a Higgs SU(2)i triplet A in the remainder of 
this section. It is straightforward to replace such a Higgs field by the product of two 
hu doublets which couple to the square of the lepton doublet L via an SU(2)i triplet 
messenger field. The type II seesaw mechanism seems particularly suited for our purposes 
as we wish to obtain the NLO correction to the neutrino mass matrix proportional to 
the coupling constant Cg, which arises from contracting LL to a 3 of T7. This can 
be naturally achieved by invoking a messenger field S'^ in the 3 representation which 
leads to the renormalizable term LLT/^. The U(l)y and U(l)/j charges required for the 
existence of this term are given in Tab. 2. Demanding a bilinear mass term fixes the 
charge assignments of S which, in turn, can couple to the Higgs field A and the flavon ^. 
A similar coupling to the other flavon is forbidden by the T^ family symmetry. We are 
therefore forced to introduce another pair of messengers f2, Qf^, which allows to generate 
the second leading order contribution to M^, proportional to ?/2, cf. Eq. (6). Their charge 
assignments are again given in Tab. 2. 

With this particle content, the renormalizable neutrino superpotential reads 

Wl^"" ~ LLS^ + S^A + SS"(Ms + + ^) 

+ LifVt + L^VL + VL^AL + nVL\Mn + V + (16) 

+ LE^fi'^ + , 

where we have suppressed all dimensionless coupling constants. As already discussed for 
the charged lepton sector, the masses of the T7 triplet messengers receive corrections from 
the flavon vevs. These will break the tribimaximal structure at higher order. Ignoring 
the corrections to the messenger masses as well as the two operators in the third line 
of Eq. (16) results in the effective superpotential of Eq. (4). The corresponding diagrams 
are sketched in Fig. 3, where the right diagram contributes to the term proportional to 
^2, while the term proportional to yi arises from both diagrams. 

We point out that the left diagram in Fig. 3 is needed only for generating the desired 
trimaximal NLO correction while suppressing other unwanted NLO terms. In other 
words, if we wanted to end up with tribimaximal neutrino mixing, the diagram with 



11 



A 



<^ A 



Figure 4: Two diagrams contributing to the NLO neutrino superpotential. The diagram on the left 
gives the desired trimaximal correction, while the diagram on the right should be suppressed in order 
to avoid large deviations of the solar angle from its tribimaximal value of d\i ~ 35.3°. 



the messenger would have been perfectly sufficient. In that case we would have to 
impose the condition (i^, ((/?) ^ Af^, similar to the situation in the charged lepton 
sector. However, since accurate tribimaximal mixing is ruled out by the Daya Bay and 
Reno measurements, we are forced to consider sizable NLO corrections. Due to the 
structure of the diagram with the Q messenger, the NLO terms obtained from attaching 
a flavon to the cross representing the mass term MqQQ'^ would give rise to corrections 
that would shift the solar angle away from its experimentally allowed region. Hence, 
the trimaximal NLO correction must originate from another diagram, namely the one 
with the S messenger, more precisely, the left diagram of Fig. 3 with the fiavon 
attached to the cross. The resulting diagram is sketched in Fig. 4, where one NLO 
diagram involving the Q messenger is also shown. The NLO diagram on the right can 
be suppressed compared to the NLO diagram on the left by assuming a hierarchy in the 
messenger masses. This can be parameterized by 

Ms ~ e^Mn , (17) 

where A; is a positive integer and e an expansion parameter around 0.2. With this 
hierarchy of messenger masses imposed we also have to assume a hierarchy in the two 
flavons vevs^ 

so that the leading order mass contributions derived from the two diagrams in Fig. 3 
are of similar size. Note that the diagram on the right involving the flavon ^ gives 
a subdominant contribution which is of the same tribimaximal structure as the one 
obtained form the left diagram. Finally, we need to guarantee that the trimaximal NLO 
correction of the left diagram in Fig. 4 is sizable enough to account for the observed 
value of reactor angle 6*13. This requires the correction to the mass of the S messenger 
originating from the flavon (/? to be of order e, i.e. 

\^e. (19) 



Ms 

Defining M ~ Mq we can summarize the requirements on the hierarchies as 

- e^^^+^M , ~ e'^+^M , Ms ~ e^M , M^^M , (20) 
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yielding neutrino mass contributions proportional to ~ ~ e at leading order, 
with a trimaximal NLO correction proportional to e^'^'^. In comparison, the right diagram 
in Fig. 4 gives a contribution proportional to e2fc+2 -^j^ic]^ is suppressed by a factor of e'^ 
relative to the desired NLO correction. 

With the above assumptions on the hierarchies it is possible to suppress all unwanted 
contractions in Eq. (12) while keeping the desired one proportional to Cg. This holds 
true even if we take into consideration possible higher order diagrams obtained from 
the two terms of the third line of Eq. (16). These terms can give rise to diagrams 
which involve one S messenger as well as one or two Q messengers. Multiplying the 
appropriate fiavon fields, one can easily find that the maximal contribution derived from 
the renormalizable term LH'^Q'^ is of order e^'^"'"^, while the one derived from Ti'^Q'^Q'^ is 
of order e^fc+s ^j^j^jg shows that these higher order corrections are also suppressed with 
respect to the trimaximal NLO correction by at least a factor of e'^. 

We conclude this section by mentioning that we have also checked these results numer- 
ically for the case of A; = 2. Integrating out the messengers we have first determined 
the effective neutrino mass matrix, from which we have calculated the mixing angles us- 
ing the Mixing Parameter Tools provided with the REAP package [54]. This numerical 
check confirms the results of this section, showing how an ultraviolet completion of the 
minimal Ty model can lead to deviations from tribimaximal mixing which are compati- 
ble with a sizable reactor angle as well as a solar angle that is close to its tribimaximal 
value. 



4. Vacuum alignment 

4.1. General discussion 

Any model of flavor that tries to explain a particular mixing pattern by means of a family 
symmetry which gets broken by fiavon vevs has to justify the imposed alignments. In 
this section we study the possibilities of achieving the vacuum configuration required by 
the minimal T7 model. Following the idea of [23], we demand that supersymmetry be 
unbroken at the scale of family symmetry breaking. This so-called F-term alignment 
mechanism introduces driving fields whose F-terms are set to zero, thus entailing the 
important F-term conditions. Just like the regular fiavons, driving fields are neutral 
under the SU(2)i x U(l)y gauge symmetry. With the five irreducible representations of 
T7 and no Zat shaping symmetry, one can introduce driving fields in at most five different 
representations. Hence, we obtain five distinct F-term conditions which we discuss in 
the following. Throughout this section we will drop all dimensionless coupling as well 
as overall CG coefficients. 

A driving field transforming in the 1 of T7 allows for two invariant terms in the 
renormalizable fiavon superpotential 

W^flav 3 Dr{M^ + ^ip) , (21) 
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which lead to the F-term condition 



+ + (^2) (^2) + (^3)(</^3) = . (22) 



Similarly, the driving fields Dy. in the representations r = 1', 1", 3, 3 of T7 lead to the 
conditions: 



dDv 
5Ds 



(^l)((^2) + ((?2)(¥^3> + (^3>(v'l) = 0, (23) 
(^l)('^3) + (^2)('^l) + ('?3)('^2) = 0, (24) 

/ (^1)2 + 2(^2X^3) 
+ L;2(^~3^2 + 2(^i)(^2)) 
V^^((^2)' + 2(^3)(^l)) 
(¥'2)(</'3)\ / (^l)(95l) + U}{lp2){(p2) + U;2(^3)(993) 

(9^3)(^l) + (^l)((^3)+^(^2)(^l)+U;2(^3)((^2) I = I I, (25) 
{fl){V2)) \(v?l)((^2) +^(¥'2)(¥'3) +U;2(993)((^i) 




(^l)2 + 2(^2)('^3) 

— 1 N--/ I + I ^((¥'3)' + 2(^l)(v^2)) 
^^3 \/,..^/ V^2(^^^^2 + 2(^3)(v,i))^ 



- mi^i) + m{vi)+^H^i){^2)+^{^2){^3) I = 1 1 . (26) 

One can easily see that the desired alignments given in Eq. (7) are inconsistent with 
Eq. (23) and Eq. (24). Eq. (25) also does not admit Eq. (7) as a solution. Turning to 
Eq. (26), one can substitute (<^i) = f^, (^2) = {^3) = and (v^i) = {(P2) = (v^s) = ^ 
to obtain 

-^M2V^ + vl + vl, + j^v^v^ = 0, (27) 
j^M2V^ + uvl + -^uj'^v^v^ = , (28) 
j^M2V^ + + j^ujv^v^ = . (29) 

Since = is not of physical interest, we can cancel off as common factor in the 
second and third equation, which can then be solved easily to obtain = ^ and 
= M2. Substituting this into the first equation, yields 2M| = 0, which in turn would 
require M2 = 0. Hence, the desired non-trivial alignment is also incompatible with 
Eq. (26). 

So far we have assumed that all terms in Eq. (25) and Eq. (26) are present. However, 
in extra dimensional models it is possible to envisage a setup in which the two fiavons 
if and If are localized on different branes [22,55]. Then, a given driving field could, in 
principle, couple exclusively to the cp fiavon, while another driving field could couple 
exclusively to (p. If both driving fields transform in the 3 of T7, it is straightforward 
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1 


1 
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2; 


-2a; 


-2|/ 


— a; — 2; 


- 2; 


u(i)« 














2 


2 


2 


2 



Table 3: The charge assignments of the hidden sector, x, y, z are positive integers smaller than N. The 
smallest possible value of -/V is 6 with {x,y,z) = (2,1,5) or equivalently (4,5,1), see discussion at the 
end of this section. 



(though tedious for the ahgnment of ip) to verify that the two flavons feature the desired 
ahgnments up to T7 transformed solutions. 

We conclude that the F-term conditions of Eq. (23)-Eq. (26) cannot give rise to the 
alignments of Eq. (7) without a mechanism that effectively sets some of the coupling 
constants to zero. 

4.2. Alignment through a hidden flavon sector 

One of the attractive features of the T7 model presented in Section 2 is that it does 
not require the introduction of an extra Z n shaping symmetry to arrive at tribimaximal 
mixing at leading order. Unfortunately, as we have found in Section 4.1, the flavons 
cannot be aligned without setting some coupling constants of the flavon superpotential 
to zero. If we do not want to resort to extra dimensions, we are forced to extend the 
minimal T7 symmetry by introducing a new symmetry. However, we will show in this 
section how this can be achieved without assigning non-trivial charges under this new 
symmetry to the flavons if and Lp, nor to the lepton and Higgs fields. 

To this end, it is necessary to introduce new flavons which do not couple to the leptons. 
In that sense, one might call them "hidden flavons" , constituting a sequestered sector of 
the model. This separation is achieved by assuming these flavons (together with some 
new driving fields) to be the only fields which transform non-trivially under a hidden 
7}^'^ symmetry. The flavons ip and (p of the lepton (or visible) sector can then be aligned 
with respect to those hidden flavons via orthogonality conditions as we discuss now. 

In the hidden sector, we introduce four flavons Xt d "^j C which are aligned by virtue 
of four driving flelds D^, D^, O^^, O^^. Their transformation properties are given in 
Tab. 3, where x, y, z are positive integers smaller than N. The minimal Z^*^ symmetry 
and the corresponding values for x, y and z will be determined at the end of this section. 
The resulting terms of the renormalizable hidden flavon superpotential read 

= 4 + XX) + D^iji; + O.xl + O • (30) 
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Setting the F-terms of the driving fields to zero gives rise to the conditions 



{X2)\ 




iXs) 


{0 + [ ^ 


{Xi)J 










-{^l){^2), 




(31) 
(32) 

(Xl)(Cl)+(X2)(C2)+(X3)(C3)=0, (33) 

(V'l) (Ci) + (V'2) (C2) + (V'3) (Ca) = . (34) 



These conditions can be solved exactly to give a total of 21 sets of solutions for the 
hidden flavon alignments. A brief description of how these can be derived is outlined 
below. 

First, we consider Eq. (32). Without loss of generality we can distinguish two cases: 
= and (^1) 7^ 0. In the former case one quickly sees that it requires vanishing 
values for all the components of (■0), thus giving only the trivial vacuum. For (-01) 7^ 0, 
we obtain three solutions, {ip) oc (1,1,1)'^, {ip) (1, cu^, cu)'^ and {tp) oc (1, cj, cj^)'^, which 
are related to each other by the T7 transformation d. We remark that a c transformation 
applied to these three solutions does not modify the alignments but only multiplies them 
with an irrelevant overall phase. 

Next, we turn to Eq. (31). Again there are two cases to distinguish, (xi) = and 
(Xi) 7^ 0- Foi" the former case the solution is easily found to be (x) oc (0, 0, 1)^. With 
(Xi) 7^ on the other hand, it is possible to show that there exist six different solutions, 
all of which are related to the previous solution by the Tj symmetry transformation c: 
(x) oc c^(0, 0, 1)"^, where A; = 0, 1, 6. One can check that applying a d transformation 
to any of these seven solutions does not yield new alignments. 

Combining our results so far, we have 21 different pairs of alignments (%) and {ijj). Being 
related by the 21 different symmetry transformations, they are physically equivalent 
and one can choose one pair without loss of generality. Having fixed the alignments 
of (x) and {ip), Eq. (33) and Eq. (34) determine the alignment of {() uniquely to be 
perpendicular to both (x) and (ip). This shows that in total there exist 21 sets of 
solutions for the hidden flavon alignments. For convenience we choose the following 
simple set 

(X) « , (V') cx , (C) oc 1^-1 j . (35) 

Let us now turn to the discussion of aligning the original flavon flelds If and (/? which 
enter in the lepton sector. Their alignment is dictated by new driving flelds which couple 
the pre-aligned hidden flavons to (p and (p. The transformation properties of the extra 
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1' 


'7' hid 


— x 


— X 


—z 


—z 


U(l)ii 
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2 



Table 4: Charge assignments of the driving fields coupling the hidden flavons x, C to the flavons ip, (fi. 



driving fields O^^, O^^ and are listed in Tab. 4. With these assignments, 

the resulting renormalizable flavon superpotential consisting of all T7 and Z^'^ invariant 
terms reads 



The derived F-term conditions take the form 

(Xl)(<?l) + (X2)(^2) + (X3)(^3) = (X3)(^3) = , 
(Xl)(<?3) + (X2)(^l) + (X3)(^2) = (X3)(^2) = , 

(Ci)(^i) + (C2) W + (C3) W = (Ci>((^i) - (</:'2)) = , 

{Cz){Vl) + (Cl)(¥'2) + (C2)(¥'3) = (Cl)((</'2) - W) = , 



yielding the desired alignments, cf. Eq. (7), 



'1' 

(^) « I I , (<^) a I 1 
.0. 



(36) 

(37) 
(38) 
(39) 
(40) 



(41) 



Finally, we need to choose an appropriate Z^"^, i.e. suitable values for A^, x, y and so 
that no additional renormalizable terms, other than the ones discussed in this section, 
occur in the flavon superpotential of Eq. (30) and Eq. (36). Through a systematic search, 
we find that the smallest Z^^ that we can have is a Zg'*^ [although Eq. (30) alone can be 
obtained using a 1}^^]. For a Zj?''^, we can have [x, y, z) as (2, 1, 5) or (4, 5, 1). Note that 
these two choices are equivalent solutions since the elements of the group corresponding 
to one set are the negatives of the elements corresponding to the other. 
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5. Conclusions 

In the present work, we discussed a minimal model of neutrino flavor. In this context, 
the term "minimal" refers to the order of our flavor group T7, and to the number of 
flavon flelds. It is interesting to note that the Z2 x Z2 Klein symmetry is not a subgroup 
of T7, but arises completely accidentally. As such, it is - to the best of our knowledge - 
the flrst indirect model in which the flavon flelds appear linearly (and not quadratically) 
in the leading order structure of the model. Another striking feature of our model is the 
absence of any Abelian shaping symmetries like Z^v or U(l). As a consequence, our flavor 
group T7 with its 21 elements is smaller than any symmetry ^4 x Zj^ that has long been 
a paradigm for model building before Daya Bay's and Reno's measurements of a large 
reactor angle ^^13 ~ 9°. Our model predicts tribimaximal mixing at leading order, and 
we achieve a sizable 613 by considering next-to-leading order corrections to the super- 
potential. It turns out that models in full agreement with experiment can be obtained 
for some generic values of the flavon vevs and coefficients in the superpotential. More 
elegantly, however, one can isolate the ncxt-to-lcading order contribution proportional 
to the coefficient C5 . Being of trimaximal type, this correction drives ^^13 to sizable val- 
ues while leaving the solar angle very close to its tribimaximal value. The atmospheric 
mixing angle is correlated to the size of ^13 and stays within the allowed 3a region. In 
Section 3, we embedded our effective theory into a renormalizable model that naturally 
suppresses all next-to-lcading order contributions except the trimaximal one which is 
responsible for the large reactor angle. In the last section, we discussed how the flavon 
flelds can be dynamically aligned to yield the required symmetry breaking pattern. The 
alignment of the flavon flelds necessitated the introduction of an Abelian shaping sym- 
metry and so-called driving flelds. However, one should note that the F-term alignment 
we considered in Section 4 is only one possible option and that the details of an elegant 
and simple mechanism of vacuum stabilization are still lurking in the dark. 
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A. Generators and Clebsch-Gordan coefficients of Tj 

In this appendix we list the relevant group theory of T7, which is sometimes also called 
the Frobenius group Z7 x Z3. The group is obtained from two generators c, d obeying 
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Table 5: The T-j generators of the five irreducible representations in the basis where d is diagonal. Here 
we have defined lo — and 77 — . 



the presentation, see e.g. [26], 

<c,d\c^ = d^ = 1, d-^cd = c^> . (42) 

It has 21 elements and five irreducible representations, namely 1, 1', 1", 3, and 3. 
A pair of triplet generators satisfying the presentation in Eq. (42) can be found e.g. 
in [14,26,27], where the order-seven generator c is diagonal. For the purpose of the 
model in the present paper, it is however more convenient to work in a basis with a 
diagonal order-three generator d. Our choice of T7 generators for all five irreducible 
representations is listed in Tab. 5, where we have defined 

a; = e 3 , rj = e ^ . (43) 

Although the order-seven generators of the triplet representations look rather involved, 
the Clebsch-Gordan coefficients take a relatively simple form. Omitting the trivial prod- 
ucts, i.e. those involving the singlet 1 as well as products of only one-dimensional ir- 
reducible representations, the product rules are reported below. We use the convention 
that the components of the first representation of any given product a (8) b are denoted 
by ttj while we use 6j for the components of the second representation. The subscripts 
s and a stand for "symmetric" and "anti-symmetric", respectively. Note that we have 
also included the normalization factors of the Clebsch-Gordan coefficients. 



1' ® 3 : 3 = ai 63 , 1" ® 3 : 3 = ai 61 , (44) 
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3 : 3 = ai 6i , 1" ® 3 : 3 = ai 63 , (45) 
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1 = 
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1" = 
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3 = 
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V3 



'1 + 0262 + 0363) , (46) 
^2 + a2h + a^bi) , (47) 

63 + 0261 + 0362) , (48) 
^aibi + ^^0262 + i^cts&aX 

0163 + a;^a26i + ^0362 , (49) 
0162 + 6^^0263 + oua^bi J 

'aihi + 6^0262 + ^'^CL-ih^ 
^ I 0361 + u;ai62 + ^^0263 I , (50) 
.0261 + u;a362 + oj'^aih^ 



aibi + 0263 + a.362 
3 ® 3 : 3, = ^ I a;(ai62 + 0261 + 0363) | , (51) 

a;^(ai63 + 0262 + ^361 



1 

V6 



^2ai6i - 0263 - 0362^ 

20262 - 0361 - 0163 I , (52) 
20363 - 0162 - 0261, 



/ 0263 - 0362 

3a = ^ 0361 - 0163 I . (53) 
\a162 - 0261^ 
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